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Abstract 

In reaction rate theory, in production-destruction type models and in reaction-diffusion problems when 
the total derivatives are replaced by fractional derivatives the solutions are obtained in terms of Mittag- 
Leffler functions and their generalizations. When fractional calculus enters into the picture the solutions 
of these problems, usually available in terms of generalized hypergeometric functions, switch to Mittag- 
Leffler functions and their generalizations into Wright functions and subsequently into Fox functions. In this 
paper, connections are established among generalized Mittag-Leffler functions, Mathai's pathway model, 
Tsallis statistics, Beck-Cohen superstatistics, and among corresponding entropic measures. The Mittag- 
Leffler function, for large values of the parameter, approaches a power-law. For values of the parameter 
close to zero, the Mittag-Leffler function behaves like a stretched exponential. The Mittag-Leffler function 
is a generalization of the exponential function and represents a deviation from the exponential paradigm 
whenever it shows up in solution of physical problems. The paper elucidates the relation between analytic 
representations of the q-exponential function that is fundamental to Tsallis statistics, Mittag-Leffler, Wright, 
and Fox functions, respectively, utilizing Mellin-Barnes integral representations. 

1. Introduction 

Fundamental laws of physics are written as equations for the time evolution of a quantity x(t), 

dx{t) 



dt 



Axit), (1.1) 



where if A is limited to a linear operator we have Maxwell's equation or Schroedinger equation, or it could 
be Newton's law of motion or Einstein's equations for geodesies if A may also be a nonlinear operator [8]. 
When A is linear then the mathematical solution is 



x{t) = xoe-^* (1.2) 

where Xq is the initial value at t = 0. Taking the simplest case of (1.1), as used in reaction rate theory 
([8], [4]), if the number density at time t of the i-th particle is Ni{t) and if the number of particles produced 
or the production rate is proportional to Ni{t) then the reaction equation is 



1 



^ = hNit), > 

deleting i for convenience. If the decay rate is also proportional to N{t) then the corresponding equation is 

dN{t) 
dt 

Then the residual effect in a production-destruction mechanism is of the form 



-k2N{t), k2>0 

fJLL 

uction-destructio 

dN{t) 



dt 



= -cN{t), c> ^ N{t) -No = -c J N{t)dt (1.3) 



if the destruction rate dominates so that the production-destruction model is a decaying model. Such 
production-destruction models abound in scientific disciplines. If the total integral or the total derivative in 
(1.3) is replaced by a fractional integral then we have 

N{t) -No = -c" o£>t""^(t), (1-4) 
where qD^'^ is the Riemann-Liouville fractional integral operator ([5], [7]) defined by 

oD^fit) = - ur-'f{u)du (1.5) 

for 5i(z/) > where 5i(-) denotes the real part of (•), c is replaced by c'' for convenience, and then the solution 
of (1.4) goes into the category of a Mittag-Lefliler function [4], namely, 

where Ei,{-) is the Mittag-Leffler function. The generalized Mittag-Leffler function is defined as 



N{t) = No J2 r{i+i) E'^i-i^r^ (1-6) 



and when r(7) is defined, it has the following Mellin-Barnes representation 

1 V- r(-, + k) 

!r(/3H 

1 1 r(s)r(7 - s) 



r(7) 2m J,_i^ r(/3 - as) 



{-zy'ds (1.9) 



for < c < 5i(7),3?(7) > 0,i = \/^. Some special cases of the generalized Mittag-Leffler function are the 
following: 

oo ^ 

El,{z) = E^,,{z) = Y: fZT-T-.^'= (1.10) 



r(/3 + ka) 

(1 + ka) 



fe=0 



K,=E.A^) = E^iz)^Y.^,^^ (1.11) 



and when a = 1 we have 

E^{z) = Ei{z) = e\ (1.12) 
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Thus the Mittag-Lcfflcr function can be looked upon as an extension of the exponential function. The 
Mellin-Barnes representation in (1.9) is a special case of the Mellin-Barnes representation of the Wright's 
function [3] which is defined as 



o'ipqiz) =p^q 



iaj,aj),j=l,...,p' 
{bi,0j),j=l,...,q _ 

Y pc-\-ioo 



OO 



(-z)-Ms 



(1.13) 
(1.14) 



where < c < niini<j<p Jfi with aj,j = l,...,p and bj,j = l,...,q being complex quantities and 

aj > 0,j — 1, ...,p and Pj > 0,j — 1, ...,q being real quantities. Observe from (1.14) that Wright's function 
is a special case of the H-function [7] and the H-function is defined as the following Mellin-Barnes integral 



(oi,ai),...,(op,ap) 
(6i,/3i),...,(6„/3,) 



1 



c+ioo 



s)z *ds 



where 



p,q — P,q 

{nu + m} {iru r(i - - ajs)} 



(1.15) 



(1.16) 



{nUm+i r(i - b, - p,s)} {uU+i + «.^)} ' 

for maxi<j<TO ''^ < c < mini<j<„ where aj,j = 1, ...,p and bj,j = 1, q are complex quantities, 

Uj > 0,j = 1, ...,p and Pj > 0,j = 1, ...,q re real quantities. Existence conditions and various contours may 
be seen from books on the H-function, for example, ([5], [7]). Equations (1.8), (1-14), and (1.15) establish 
the relations between Mittag-LefFler, Wright, and Fox functions to be used in the following for q-exponential 
function in Tsallis statistics, Beck-Cohen superstatistics, and their link to Mathai's pathway model. 

1.1. Extension of the reaction rate model and Mittag-Leffler function 

The reaction rate model in (1.3) can be extended in various directions ([8], [4]). For example, if No is 
replaced by Nof{t) where f{t) is a general integrable function on the finite interval [0,6] then it is easy to 
see that for the solution of the equation 



Nit) - Nofit) = -c" oD^Nit) 



there holds the formula 



N{t) = cNo f Hl'l 
Jo 





(-i.i) 








(-i,l),(0,i')_ 



(1.17) 



(1.18) 



Some special cases are the following: Let i/ > 0, p > 0, c > 0. Then for the solution of the fractional equation 

N{t) - NotP-^ = -c" oD-^Nit) 



there holds the formula 

N{t) = NoTipy-^E^^pi-icty). 
Let c > 0, 1/ > 0, /U > 0. Then for the solution of the fractional equation 



(1.19) 
(1.20) 



N{t) - Not^-'E2J-{ctr] = -d' oDrN{t) (1.21) 

there holds the formula 

N{t) = Not''-'E2+'[-{ctn (1-22) 

1.2. Fractional partial differential equations and Mittag-Leffler functions 

In a series of papers, see for example, ([2], [4]) it is illustrated that the solutions of certain fractional 
partial differential equations, resulting from fractional diffusion problems, are available in terms of Mittag- 
LefHer functions. For example, consider the equation 

oD'^Nix, t) - Y^^^^i^) = -^'^^(^' ^) (1-23) 



with initial conditions 



oD'^-''N{x,t)\t=o = 0,k = l,. 



..,n 



where n = [J?(z^)] + 1, c'' is the diffusion constant, 6{x) is the Dirac's delta function and [5R(i')] is the integer 
part of 5ft(i^). The solution of (1.23), by taking Laplace transform with respect to t nd Fourier transform 
with respect to x and then inverting, can be shown to be of the form 



N{x,t) = r^f 2 

which in special cases reduce to Mittag-Leffler functions. 





(i-f,'') " 




(0,1), 



(1.24) 



This paper is organized as follows: Section 2 gives the classical special function technique of eliminating 
the upper or lower parameters from a general hypergeometric series. Section 3 establishes the pathways of 
going from Mittag-Leffler function to pathway models to Tsallis statistics and Beck-Cohen superstatistics 
through the parameter elimination technique. Section 4 gives representations of Mittag-Leffler functions and 
pathway model in terms of H-functions and then gives a pathway to go from a Mittag-Leffler function to the 
pathway model through H-function. 

2. A classical special function technique 

The well established technique of eliminating a numerator or denominator parameter from a general 
hypergeometric function in the classical theory of special functions, is the following: For convenience, we 
will illustrate it on a confluent hypergeometric function. 



oo / \ ^ 

iFi(a; 6; ^) = ^ (a)„ = a{a + l)...(a + m - 1), a 7^ 0, (a)o = 1. (2.1) 

fe=0 ^ ' 

Observe that 



{a)k _a{a + l) {a + m-l) 
a'' a a a 

= 1(1 + ^)(1 + ^)...(1 + ^^^^) ^ 1 as a ^ 00 (2.2) 

for all finite k. Similarly 

— > 1 when b 00. (2-3) 

Therefore 
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lim iFi{a;b;-) = oFi{ ■,b;z) (2.4) 

a^oo a 

which is a Bcsscl function. Thus we can go from a confluent hypergeometric function to a Bessel function 
through this process. Similarly 

lim iFi(a;6;62) = iFo(a; : z) = {1 - z)-", \z\ < 1. (2.5) 

6— »oo 

Thus, from a confluent hypergeometric function we can go to a binomial function. Further, 

lim oFi{ ; b; bz) = oi=o( ; ■,z) = e' (2.6) 

6— »oo 

and 

lim iFo(a; ■,-)=oFo{ ; ;-^) = e^ (2.7) 
o— >oo a 

Thus, we can go from a Bessel function as well as from a binomial function to an exponential function. These 
two results can be stated in a slightly different form as follows: 

limoi^i( : ^)=e-" (2.8) 

q^i ^ q-1 q-1' 

and 

limii=^o(^-; ;-{q-l)z) = \im[l + {q-l)z]-^ =e-'. (2.9) 

q^l q — I 9^1 

Equation (2.9) is the starting point of Tsallis statistics [9]. The left side in (2.9) is the g-exponential function 
of Tsallis, namely, 

[l + {q-l)zr^ =exp^(-z) 

= e"^ when q ^ 1. (2.10) 

We consider a more general form of (2.9) given by 

lim c\x\^iFo{^^; ; ~a{q - l)\xf) = lim c\xp[l + a{q - l)|x| 

= c|a;Pe-"''l^l',a > 0,r? > 0, -oo < a; < oo. (2.11) 
Thus, we obtain Mathai's pathway model [6], namely, 

f{x) =c\xy[l + a{q-l)\x\^r^ ,a>{),'n>0 (2.12) 

where c is the normalizing constant. If = 1,7 = 0,a = 1,(5 = l,x > then (2.12) gives Tsallis statistics 
[9]. For g > l,a = l,ry = l,a; > in (2.12) we get Beck-Cohen superstatistics [1]. 

It may be observed that (2.7), which is the binomial form or iFq going to exponential, is exploited to 
produce the pathway model, Tsallis statistics and Beck-Cohen superstatistics whereas (2.6), which is the 
Bessel function form going to exponential is not yet exploited. This produces a rich variety of applicable 
functions as will be shown in a forthcoming paper. 

3. Connection of Mittag-Leffler function to the pathway model 

In order to see the connection, let us recall the Mellin-Barnes representation of the generalized Mittag- 
Leffler function. 
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for 3^(7) > 0,i = 3fi(/3) > 0. Let us examine the situation when |/?| 00. We will state two basic 

results as lemmas. 

Lemma 3.1 When — > 00 and a > is finite then 



Proof. The Stirling's formula is given by 

r(z + a) « v^^^+"-^e-^ for \z\ 00 (3.3) 
and a is a bounded quantity. Now, applying Stirling's formula wc have 



2 c 



-0 



p'^'i-bzY'P""" = {-bz^)~"- (3.4) 



Lemma 3.2. For 3?(7) > 0, 3t{(3) > 0, 



1 1 pc+ioo 



/c=0 



(3.5) 
(3.6) 



Proof. After applying Lemma 3.1 we can evaluate the contour integral in (3.5) by using the residue theorem 
at the poles of T{s) to obtain the series form in (3.6). For \bz^\ > 1 we obtain the series form as analytic 
continuation or by evaluating the integral in (3.5) as the sum of the residues at the poles of r(7 — s). This 
will be equal to the following: 

lim VmEl^{-bz') = {hz')-'^[l + {hz')-']-'^ (3.7) 

for \bz^\ > 1 which will be of the same form as in (3.6). Replacing b by a{q — 1) and 7 by ^^j- we have the 
following: 



lim c\zrmE2^l;'-'\-a{q-l)\z\') 

= c\z\^[l+a{q-l)\zf]-^. (3.8) 

Observe that (3.8) is nothing but Mathai's pathway model [6] for a > 0, (5 > 0, 77 > with c being the 
normalizing constant, for both q > 1 and q < 1. Then when q 1 (3.8) will reduce to the exponential form. 
That is, 

lim^^^lim c\z\^T{/3)El(^^''-^^[-a{q-l)\zf]=c\z\^e-'''\'\\ (3.9) 
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Note from (3.8) that when q > I the functional form in (3.8) remains the same whether (5 > or (5 < with 
— oo < z < oo. When g < 1 the support in (3.8) will be different for 5 > and (5 < if f{x) is to remain as 
a statistical density. 

For 7 = 0,o=l,(5=l,?7 = l,^;>0we have Tsallis statistics [9] coming from (3.8) for both the cases 
q > 1 and q < 1. For q>l,a=l,5=l,rj^l,z>0 we have Beck-Cohen superstatistics [1] coming from 
(3.8). Observe that the limiting process from (3.5) to (3.6) holds for 7 = 1 or 7 = 1 and a = 1 also. Thus 
the special cases of Mittag-LefHer function are also covered. 

Thus, through (3 a pathway is created to go from a generalized Mittag-LefHer function to Mathai's 
pathway model and then to Tsallis statistics and Beck-Cohen superstatistics. If f3 is real then as f3 becomes 
larger and larger then 



c\z\ 



y(.-i)(_«(g_i)l,l 



goes closer and closer to the pathway model in (3.8). In other words, a pathway is created through /3 to 
go from a Mittag-LefHer function to the pathway model to Tsallis statistics and Beck-Cohen superstatistics. 
Thus for large real value of (3 or for large value of |/3|, 



c|^rr(/3)£;^f^«-') [-a(g - 1)|^|*] « c|^r [1 + a(g - l)\zf]-^. 
4. Connections through the H-function 



(3.10) 



Recalling the generalized Mittag-Leffler function from (3.1) and representing it in terms of a H-function 
we have the following: 



J c—ioo -^(/^ 



(7 - s) I 



as) 



-(z/3*)*]-Ms 



r(7) 



^1,2 





(1-7,1) 


-{z(3f)' 






(0,l),(l-/3,a)_ 



(4.1) 

(4.2) 



From the limiting process discussed in (3.2) we have the following result: 
Lemma 4.1 For 3?(/3) > 0, ^{j) > 0, 



|/3Hoor(7) 



51(1-7,1) 



rrl.l 

-^1,1 



(0,l),(l-/3,a) 

5|(1-7,1)' 
1(0,1) 



r(7) 

1 -1 /■C-l-lOO 



/■c+ioo 

r(7) 2^ 
[1 - /]-^. 



Ms 



(4.3) 

(4.4) 
(4.5) 



Therefore we have the following theorem: 

Theorem 4.1 For 3?(/3) > 0, > 0, x > 0, a > 0, g > 1, c> 



lim c^^x-'<(-^)[-a(,-l)/3fx]^ 
= cx^[l + a{q - l)x^]~^ . 



(4.6) 



Observe that the right side in (4.6) is the pathway model [6] for a; > from where one has TsalUs statistics, 
Beck-Cohen superstatistics and power law where the constant c can act as the normalizing constant to create 
a statistical density in the right side of (4.6). 

We will write the right side in (4.6) as a H-function and then establish a connection between Mittag- 
LefHer function and the pathway model through the H-function. To this end, the practical procedure is 
to look at the Mellin transform of the right side of (4.6) and then write it as a H-function by taking the 
inverse Mellin transform. The Mellin transform of the right side of (4.6), denoted by Mf{s), is given by the 
following: 

POO 

Mf{s)= / cx'^+^-'^[l + a{q-l)x^]--^dx,a>0,q>l,ri>0,S>0 
Jo 



_n 7+sJ 

-1 5 



(4.7) 



for 5J ("^"T^j > 0, — ^^7^) > 0. Note that if c is the normalizing constant for the density f{x) then 
by putting s = 1 the right side of (4.7) must be 1. Therefore, 



Mf{s) 



7+1 
S 



7 + s 



q-l 



i(g-l)]-t. 



Hence the right side of (4.6) is available as the inverse Mellin transform of (4.9). That is, 



(4.9) 



where 



f{x) =cx^[l + a{q-l)x^\-^, 

for a > 0, > 0, g > 1, ^ > 0, a; > 



7 -I- s 



7 + s 



q-l 



[a{q-\)]-s 



H. 



1,1 



x{a{q-l))^ 



= lim c 



_ J[a(g-l)]^r(^) 

'"r(2+i)r(^-2+i) 

for g > 1, 5 > 0, 77 > 0, a > 0, 3?(7 1) > 0, - 



xia{q-l))^ 



Ms 



7± 
S 



> 0. 



(4.10) 



(4.11) 
(4.12) 

(4.13) 
(4.14) 



Remark 4.1. From (4.13) and (4.10) it can be noted that as the parameter f} becomes larger and larger the 
Mittag-LefHer function in (4.13) comes closer and closer to the pathway model and eventually in the limiting 
situation both become identical. In a physical situation, if (4.10) represents the stable situation then the 
unstable neighborhoods are given by (4.13). When q ^ 1 then (4.10) goes to the exponential form. Thus 
if the exponential form is the stable situation then the pathway model of (4.1) itself models the unstable 
neighborhoods. This unstable neighborhood is farther extended by the Mittag-LefHer form in (4.13). 
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Remark 4.2. Observe that the funetional form in (4.10) remains the same whether g' > 1 or g' < 1. But 
for g < 1 or when g — > 1 the normahzing eonstant c will be different. 

Remark 4.3. When dealing with problems sneh as reaction-diffusion situations or a general produetion- 
destruction model one goes to fractional differential equations to get a better picture of the solution. Then 
we usually end up in Mittag-Leffler functions and their generalizations into Wright's function or Fox function. 
Comparison of (4.13) and (4.10) reveals that as (3 gets larger and larger the effect of fractional derivative 
becomes less and less and finally when |/3| oo the effect of taking fractional derivatives, instead of total 
derivatives, gets nullified. 
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